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ABSTRACT 

One of the puzzles associated with tidal disruption event candidates (TDEs) is that 
there is a dichotomy between the color temperatures of few x 10^ K for TDEs discovered 
with optical and UV telescopes, and the color temperatures of few x 10^ — 10® K for 
TDEs discovered with X-ray satellites. Here we propose that high-temperature TDEs 
are produced when the tidal debris of a disrupted star self-intersects relatively close to 
the supermassive black hole, in contrast to the more distant self-intersection that leads 
to lower color temperatures. In particular, we note from simple ballistic considerations 
that greater apsidal precession in an orbit is the key to closer self-intersection. Thus 
larger values of j3, the ratio of the tidal radius to the pericenter distance of the initial 
orbit, are more likely to lead to higher temperatures of more compact disks which are 
super-Eddington and geometrically and optically thick. Eor a given star and /3, apsidal 
precession also increases for larger black hole masses, but larger black hole masses imply 
a lower temperature at the Eddington luminosity. Thus the expected dependence of the 
temperature on the mass of the black hole is non-monotonic. We find that in order 
to produce a soft X-ray temperature TDE, a deep plunging stellar orbit with /? > 3 
is needed and a black hole mass of < 5 x 10®Mq is favored. Although observations 
of TDEs are comparatively scarce and are likely dominated by selection effects, it is 
encouraging that both expectations are consistent with current data. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: nuclei — 
relativistic processes — stars: kinematics and dynamics — X-rays: bursts 



1. Introduction 


The discovery of roughly two dozen TDEs with X-ray and optical/UV telescopes (e.g., [Bade 


et al.||1996t Komossa &: Greiner||1999; Esquej et al.||2008 van Velzen et al.||2011 Gezari et al.||2009 

2012 

Genko et al. 

2012 

Saxton et al. 20121 Maksym et al. 

2013 

Ghornock et al. 2014 

Arcavi 


et al. 2014; Holoien et al. 2014, 2015) has afforded us the opportunity to study accretion over 


a wide range of rates and is also promising for the discovery and characterization of otherwise 
quiescent supermassive black holes (SMBHs) ( Loeb &; Ulmer||1997 ). Standard treatments predict a 
temperature that is few x 10^ K near the peak of the burst and decreases with the luminosity (Rees 


1988; Phinney 1989; Cannizzo et al. 1990 Ulmer 1999). However, although X-ray detected TDEs 


reach this temperature, the TDEs discovered via optical/UV observations have lower temperatures 
few X 10^ K, and those temperatures remain nearly steady even as the luminosity of the sources drop 
by 1 — 2 orders of magnitude (e.g., Gezari et al.||2012 Arcavi et al.||2014 ). This has been explained 
as a consequence of either an optically thick shroud of gas at many times the radius of the disk 
(Loeb & Ulmer 1997) or a wind from the disk ( Strubbe &: Quataert] 2009 Millerj 2015 Metzger & 


Stone 2015). Why do the processes which reduce the temperature in optical/UV detected TDEs 
fail to operate for the X-ray detected TDEs that have temperatures of few x 10^ — 10® K? 

The ranges in length and time scales required for full hydrodynamic simulations of debris 
circularization in tidal disruptions of a main-sequence star by a SMBH from a marginally bound 
orbit mean that such simulations are currently computationally infeasible. The main reason is that 
such simulations need to be able to follow the debris orbit with radius r ~ 100 — lOOOi?^, and at 
the same time need to resolve the thickness of the stream which, if is self-gravitating, is much less 
than r (Kochanek 1994; Guillochon et al(]|2014 ). The addition of realistic cooling mechanisms and 
general relativistic hydrodynamics would make simulations even more time-consuming. However, 


studies of stars disrupted by SMBHs from initially bound orbits (Bonnerot et al. 2015 Hayasaki 


et al.|20T5 ) and white dwarfs disrupted by intermediate mass black holes (Ramirez-Ruiz Sz Rosswog 
2009; Shiokawa et al.|2015 ) have revealed that debris circularization is more difficult than previously 
thought. The nozzle shock and the instabilities at pericenter are not strong enough to circularize 
the debris quickly (Kochanek |1994 Guillochon et al. 2014 Shiokawa et al. 2015). Instead, tidal 
stream intersection is the most effective way to produce shocks and dissipate the debris orbital 


energy (Shiokawa et al. 2015; Bonnerot et al. 2015 Hayasaki et al. 2015). This intersection is a 


result of the orbital apsidal precession of the debris on the same plane around the SMBH. However, 
if apsidal precession is small and intersection happens far from the pericenter of the initial orbit. 


then the collision is mild and circularization can be delayed substantially (see, e.g., Shiokawa et al. 


2015 and the Newtonian simulation of Bonnerot et al. 2015). In this case, most of the matter 


would take much longer to accrete onto the black hole than it would in the standard picture. As a 
consequence, the luminosity and the disk temperature would both be less than those seen in X-ray 
detected TDEs. 


Here we propose that the high temperatures of some TDEs can be explained in a picture in 
which the debris disk is small due to relatively large apsidal precession. In particular, we argue 
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that if the disk is small enough that the initial temperature exceeds ~ 10^ K, then (unlike when 
T = few X 10^ K) the opacity is only weakly sensitive to the temperature, which eliminates the 
strong dependence of wind rates on temperatures that features in one explanation for the lower 
temperatures seen in optical/UV TDEs (Miller 2015). In Section 2 we perform calculations of 
ballistic motion to explore how debris stream intersection and the consequent disk size depend on 
the black hole mass M and the stellar orbit penetration parameter j3. In Section 3 we calculate 
the disk temperature, where we show that high-/? encounters, particularly with black holes of 
M ~ 10® Mq, are the best candidates for TDEs with temperatures > few x 10® K. Eurther 
discussions and some observational considerations are given in Section 4. 


2. Tidal Stream Intersection: First Order Calcnlations 


As the nozzle shock is weak, the trajectory of the debris is nearly ballistic until the tidal 
stream self-intersects. Although some self-intersection would occur even in Newtonian gravity 
because there is a spread in debris binding energy, the dynamics of the intersection are dominated 
in our case by general relativistic pericenter precession. If we assume for simplicity a non-rotating 
(Schwarzschild) black hole (rotational corrections are mild unless the pericenter is very close to the 
hole), then the precession angle (j) over a single orbit is, to first order, 

(p = 67rM/(o(l — e^)), (1) 


(Misner et al. 1973), where a is the semi-major axis of the orbit, and e is the orbital eccentricity. 
Here and henceforth we use geometrized units in which G = c = 1. If the orbit is close to the hole, 
the exact precession rate is larger than the rate given by this expression, but the differences are 
mild for a(l — e) > lOM. In reality the ellipse precesses continuously, but for the high-eccentricity 
orbits we consider almost all the precession occurs at pericenter. We therefore treat the debris orbit 
as a closed ellipse with an instantaneous shift of (p at pericenter passage, as illustrated in Eigure[TJ 


The head of the stream essentially follows the most bound debris orbit, and after passing 
pericenter it intersects the trailing part of the incoming stream (the shifted orbit is represented 
by the blue ellipse in Eigure[^. The incoming stream lies between the most bound debris orbit 
(the black ellipse) and the marginally bound debris orbit (the black dotted curve). If the initial 
trajectory of the star is nearly parabolic, intersection typically happens in a time after disruption 


of I - I.5Pmb ( 

Shiokawa et al. 

2015 

Hayasaki et al. 

2015 

Bonnerot et al. 

2015 

Guillochon & 

Ramirez-Ruiz 

2015 

), where Pmb is the orbital period of the most bound orbit. 


We now study the dynamics of the most bound orbit for the disruption of a star of mass 
m* = 10® and a radius = 10® R^aRq, where Rq is the radius of the Sun, by a SMBH 

with mass M = 10® MqMq. 
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2.1. The most bound orbit 

A star on an initially parabolic orbit has zero specific binding energy: 

E = 1/2 - M/Rt = 0, (2) 

where Rt ~ R^{M is the tidal radius and vt is the orbital speed at Rt- If we neglect the 
rotation of the star and the redistribution of energy during the compression / rebound processes, 
the most bound orbit has energy: 


-E'mb = 1/2 v'^ — M/{Rt — R*) ~ —MR^/R^, 


( 3 ) 


because Rt S> R*- The semi-major axis and eccentricity of this orbit are 

Omb = Rt/‘^R* ~ 3.5 X 10^“^ cm 


( 4 ) 


and 

emb = 1 - Rp/amh 1 - 2//3 X = 1 - 0.02m^|'Q^Mg"^^^/3“^ (5) 

Here Rp is the pericenter distance, and we have defined the penetration factor /3 = Rx/Rp- The 
orbital period of the most bound debris is: 


Pmb = 27 := 0.11 yr rI^q 


( 6 ) 


(note that we neglect the weak increase in Pmb with increasing /3 found by Guillochon & Ramirez- 

H disruption is hig 
1 year (Mg = 100). 


Ruiz 2013). Therefore, the most bound orbit in a main-sequence star-SMBH disruption is highly 

1 month (Mg = 1) and 


eccentric and its orbital period ranges between 


2.2. The intersection radius and collision angle 

As discussed previously, the narrowness of the tidal stream means that the intersection radius 
of the tidal streams is given by apsidal precession. Treating the most bound orbit as an ellipse with 
an instantaneous pericenter shift cj), a geometrical calculation shows that the original ellipse and 
the shifted ellipse intersect at a radius 

/? — (1 + emh)RT 

'“/3(l-e^bCos(<^/2))' ^ ^ 

Figure l^a) shows that this stream intersection radius changes sensitively with M and /3 (hereafter 
we use a solar-type star in calculations). For more massive holes, the tidal radius relative to the 
size of the hole decreases, so stronger apsidal precession reduces the intersection radius in units of 
the black hole gravitational radius Rg = GMjc?. For higher /3, Rp is smaller, which again implies 
stronger precession and closer intersection. 
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The intersection angle 0 of the outgoing most bound orbit with the incoming stream is given 
by 

cos 0 = l-2co«(0/2)e...b + cos,»4c (8, 

1 - 2cos((/)/2)emb + e^b 

0 is plotted as a function of M and /3 in Figure j^b). Streams usually collide with 0 ~ 40 — 160° 
if 2 > Rg. The larger the intersection angle, the more effectively the collision reduces the orbital 
energy of the debris. 


2.3. Energy dissipation at stream intersection 


Some fraction of orbital energy is lost during stream collision due to shocks and instabilities. 
The exact fraction depends on the collision angle and velocity/density contrast of the streams 
(jKochanA 1994). Here we adopt an inelastic collision model in which the outgoing and incoming 


streams have similar mass (cf. Shiokawa et al. 2015). From momentum conservation, the post¬ 


collision speed of the stream at the intersection point is 

Vf = Vi cos(0/2) , 

where the speed Vi of the streams just before collision is given by 


GM 

2Umb 


GM 1 


Ri 


+ 2"* 


(9) 


( 10 ) 


The specific energy loss in the collision is AE = ^v’f sin^(0/2). The debris forms an elliptical 
disk after collision, and the semi-major axis of this elliptical disk is 

Rj 1 


®disk — 


2 sin2(0/2) 1 + ^ cot2(0/2)' 
For deep encounters, Omb Ri and thus 

Ri 


O-disk 


2sin2(0/2)' 


( 11 ) 


( 12 ) 


We plot this characteristic size of debris disk in Figure The disk size shrinks rapidly as M or /3 
increases. This is because the energy loss is large when the intersection is close and the collision 
speed is fast. When (3 is large, the disk radius is comparable to the classical circularization radius 
Rc = 2/13 X Rt, for which the debris materials are fully circularized. 


3. Debris Disk Accretion and Temperature 

The evolution of the post-collision elliptical disk depends on further stream-stream and stream- 
disk interactions. The circularization timescale has traditionally been assumed to be several dy¬ 
namical timescales of the most bound orbit: 

^circ — n X Tjjjb, (13) 
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with n = 2 — 10 (Evans &: Kochanek 1989; Ulmer|1999 ). However, Shiokawa et al. (2015) found that 
for a white dwarf disrupted by an intermediate-mass black hole the disk is only partially circularized 
after n > 10 if the stream collision is mild. Though it is unclear how this result extrapolates to a 
higher mass ratio TDE, it is plausible that the circularization timescale will be shorter for higher /3 
where the collision is stronger. Eor our calculation we assume at Tcirc = 5Pmb the debris materials 
are largely circularized and the disk becomes stable, which corresponds to the peak of the disk 
luminosity, but our results are not qualitatively changed unless n is at least an order of magnitude 
larger than we assume. 


Erom the classical a-disk model ( Shakura &: Sunyaev||1973 ), the inflow timescale of the disk is 

(14) 


_l (H 

Tinflow- - X ( - 


-2 


X Pn 


where a is the viscosity parameter of the disk, H/R is the disk aspect ratio at a radius R, and 
Pa is the orbital period at the outer edge of a disk of size Odisk- Therefore, in deep plunging 
TDEs, because Odisk Omb; Pa Pmh, and the inflow timescale is likely to be shorter than the 
circularization / fallback timescale. In this situation the onset of accretion is therefore rapid; it can 
easily occur before the disk circularizes fully. 

Eor a TDE around a black hole with M < few x IO^Mq, the mass fallback rate is greater 
than the Eddington accretion rate for the first weeks to years. When the circularization is efficient 
due to the strong stream collisions as discussed above, the gas supply rate to the disk is also super- 
Eddington. Modern simulations of super-Eddington accretion flow (e.g., Sqdowski et al.|2015 ) show 
that such disks are geometrically and optically thick, which together with strong outflows completely 
obscure the inner region of the disk. Studies of the outflow structure and rate of super-Eddington 
accretion are far from complete, so for simplicity we assume that ~ 10% of the fallback mass flows 
out in a wind, launched from near the innermost stable circular orbit (ISCO) at the orbital speed 
there (for motivation see Blandford Sz Payne|1982 Tchekhovskoy et al.j2012l . At the circularization 
time, the wind from the power-law decay part of the fallback (starting at ~ 1.5Pmb) has reached 
out to Tout = (Tcirc “ l-5Pmb) X vk- Taking the time-averaged M^ind - 0.1 x Mfaiiback, the 
optical depth of the wind beyond ajisk is: 


r = 


r° / X , r ^wind J KesMfallback , 

/ Kesp{r)dr = / Kes-. - 3 - 

where Kes = 0.34 cm^ g“^ is the electron scattering opacity, and vk = c/2 at the Schwarzschild 


(15) 


1 

Q-disk 


1 , 
Cout 


ISCO. Using Mfaiiback(t) from Evans & Kochanek (1989), and given rout ^disk, we find that r ~ 
KesTlfaiiback/(407ruxadisk) < 1 except in the case of the most extreme super-Eddington accretion, 
for which r ~ 1. Thus the wind beyond the disk does not significantly change the observed 
temperature. 

We therefore compute the effective temperature Tefj of the system using a photospheric size 
~ fldisk: 

L = 47ra < TEdd- (16) 























Here L is the luminosity of the disk, a is the Stefan-Boltzmann constant, and we cap the luminosity 


at Eddington; this is consistent with some simulations (McKinney et al. 

2014), 

although other 

simulations find that the luminosity could be somewhat higher (Jiang et al. 

2014; 

S§,dowski et al. 


2015), which would imply higher temperatures. The accretion efficiency r] is always taken to be 


~ 0.1, which is comparable to the specific binding energy at the ISCO of a moderately spinning 
black hole. This calculation of the photosphere size is applicable to super-Eddington accretion. If 
the fallback rate at 5Pmb is sub-Eddington, we use the fallback rate to calculate the peak accretion 
power, but Teg oc is insensitive to Macc- 

We plot Tgfj in Figure]^ increases with increasing (5. To reach the T > few x 10® K seen 
in X-ray selected TDEs, it is necessary that P > 3 and M < 5 x 10 ®Mq. Note that disk radiative 
transfer effects can alter the spectrum so that the best-fit Planck temperature can be up to ~ 2 
times larger than the effective temperature (see, Table 1 of Davis et al.|2()05 ). Therefore, very large 
P may not be required to explain TDEs with T > 5 X 10® K. 


4. Discussion and Summary 


We consider the tidal disruption of a main-sequence star by a non-spinning SMBH, and in 
particular the tidal stream intersection and circularization driven by apsidal precession of the 
debris orbit. We show that because apsidal precession is greater for more massive black holes and 
deeper orbital penetration, such encounters also lead to stream intersection closer to the hole. This 
strengthens shocks, which enhances energy losses and decreases the circularization timescale. The 
resulting super-Eddington accretion disk is small and has a short inflow timescale, so the initial 
disk temperature is high. 


For a fixed black hole mass, stars disrupted in closer orbits produce higher temperature events. 
However, for hxed /3 this disk temperature does not decrease monotonically with increasing black 
hole mass; it appears that for X-ray detected TDEs it is necessary that j3 > 3 and M < 5 x 10® M©. 
This is consistent with the mass distribution of the small number of X-ray observed TDE flares 
reported in [Stone &; Metzger (2014). Such high (3 encounters are most likely to be produced when 
the orbital phase space of the stars is not depleted by disruptions. In this case, the probability of 
an event with /3 > /3o scales as l//3o- If, in contrast, the inner galactic zone is depleted, diffusion 
processes tend to produce low-/3 encounters (e.g., Alexander]2005 MacLeod et al.j2013 ). Therefore 
we expect that < 1/3 of the fully disruptive events reach a temperature high enough to produce 
substantial X-rays. Reprocessing from, e.g., dust could reduce this fraction. Nonetheless, this ratio 
is consistent with the ratio of the observed X-ray TDE rate (Donley et al. 2002) and optical/UV 
TDE rate (van Velzen Sz Farrar 2014), though both estimates bear uncertainties due to the small 


sample sizes and modeling assumptions. We note that not many TDEs have been observed at their 
peaks with X-ray telescopes, so a more complete search on X-ray TDEs will greatly improve our 
understanding. 




























There are other models proposed to explain X-ray TDEs. For example, inferred temperature 
differences among TDEs may be partially due to viewing angle dependent obscuration by the disk 


(Watarai et ah 2005; Coughlin Sz Begelman 2014) or wind (McKinney et ah 2014; S§,dowski & 


Narayan]|2015 ), leading in the extreme case to a jetted TDE when viewed down the jet (McKinney 


et al.|[2M5). If the disk is geometrically thin, unobscured, and the accretion rate can stay relatively 


constant with radius, then the disk temperature can be dominated by the temperature near disk 
inner regions (Guillochon et al. 2014). Non-thermal processes such as electron scattering and 
Comptonization can drive part of the TDE disk to have a high temperature (Li et al. 2002). 
Thermal or non-thermal X-ray emission could also arise from either the disk corona or jet, where a 
similar ambiguity exists for the origin of SgrA* radio and X-ray emission ( Falcke Sz Markoff||2000 ). 

Around a spinning black hole, the nodal precession of the debris stream can delay substantially 


the time of first stream intersection (Dai et al. 2013; Hayasaki et al. 2015; Guillochon &z Ramirez- 


Ruiz 1(2^5), and this likely also delays the onset of accretion. As debris streams may only partially 


collide, circularization can take longer and the peak luminosity will be lower than it would be 
around a non-spinning hole. However, if the initial stellar orbit is aligned with the black hole spin, 
the stream intersection will be very similar to that discussed in this paper. The extra precession 
produced by black hole spin on a retrograde orbit can draw the intersection closer to the hole. 
This produces an even smaller disk and thus promotes faster onset of accretion and increases disk 
temperature. 
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Fig. 1.— The intersection of tidal streams for different setups. The black hole is at the origin. The 
black solid ellipse is the most bound orbit in a Newtonian potential. The black dotted curve is the 
parabolic trajectory of the center of the star. The blue dashed ellipse is the ballistic trajectory of 
the most bound orbit with a periastron precession of angle cj) (Equation (1)). Panels (a), (b), and 
(c) represent the scenarios for a SMBH and a solar type star. Panel (d) represents an IMBH-white 
dwarf disruption scenario as in Shiokawa et al. (2015). Stream-stream intersection occurs where 
the black solid line and blue dashed line intersect. In a deep plunge (panel (b)) or when the black 
hole is more massive (panel (c)) the stream intersects closer to the hole. 
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(a) Stream Intersection Radius 



(b) Stream Collision Angle 



Fig. 2.— The intersection radius and collision angle of the most bound orbit with itself after apsidal 
precession. Different colors and line styles represent different /?: 1 - black solid, 2 - blue dotted, 4 
- green dashed, 8 - orange dot-dashed , and 16 - magenta long-dashed. Panel (a) shows that the 
interaction radius decreases with larger /3. Panel (b) shows that the collision angle in most cases is 
between 40 and 160°. The energy loss is more efficient when the collision angle is large. 
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Fig. 3.— The characteristic sizes of the elliptical disk after stream intersection are represented 
by the thick curves. The color and line style scheme are the same as in the previous figure. For 
comparison, the thin lines show the corresponding classical circularization radius for parabolic 
TDEs. This figure shows that the disk size greatly decreases with increasing /3. 
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Fig. 4.— The effective temperature of the debris disk. The color and line style scheme are the 
same as in Fig. 2. The light blue thin line is the temperature of a thin disk at its ISCO. This figure 
shows that soft X-ray temperature TDEs are produced by stars in deep plunging orbits with /3 > 3 
(green-dashed line and lines with higher Teg). 










